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We present a detailed investigation concerning the exciton dynamics in GaN epilayers grown on 
c-plane sapphire substrates, focussing on the exciton formation and the transition from the non- 
thermal to the thermal regime. The time-resolved kinetics of LO-phonon replicas is used to address 
the energy relaxation in the excitonic band. From ps time-resolved spectra we bring evidence for a 
long lasting non-thermal excitonic distribution which accounts for the first 50 ps. Such a behavior is 
confirmed in different experimental conditions, both when non-resonant and resonant excitation are 
used. At low excitation power density the exciton formation and their subsequent thermalization is 
dominated by impurity scattering rather than by acoustic phonon scattering. The estimate of the 
average energy of the excitons as a function of delay after the excitation pulse provides information 
on the relaxation time, which describes the evolution of the exciton population to the thermal 
regime. 



PACS numbers: 78.30.Fs, 78.55. Cr,78.47.jd,71.35.Cc 

I. INTRODUCTION 

In direct band gap semiconductors the presence of ex- 
citonic resonances in optical spectra is considered as a 
clear marker to assess the sample quality. In the paradig- 
matic GaAs based systems, excitons are usually probed 
by means of photoluminescence (PL) experiments and 
their dynamics is investigated in a restricted wave vector 
(k) interval near the T point (close to k = 0) as a conse- 
quence of the radiative recombination process which re- 
quires, in the photon emission, both energy and momen- 
tum conservation. Therefore, when an excess energy is 
used to excite the sample, the risetime of the PL signal at 
the exciton energy includes all the information concern- 
ing exciton formation and their relaxation towards the 
k«0 radiative states, making difficult to extract quan- 
titative information on each process. Even when reso- 
nant excitation is used the information on the k^O exci- 
tons can only be obtained by a careful modeling of the 
exciton dynamics.^ For semiconductors where excitons 
strongly interact with the lattice, phonon replicas (PR) 
of the k = excitonic band ( the so called zero-phonon 
line (ZPL)) are present in the PL spectra. Two possi- 
ble situations can occur depending on the exciton local- 
ization/delocalization process in the material. In disor- 
dered/defective materials - where most of the recombi- 
nations come from extrinsic contributions - excitons are 
strongly localized. The Huang-Rhys (H-R) model is then 
commonly used to describe the replica characteristics and 
to extract the H-R coefficient, which quantifies the exci- 



ton interaction with the lattice. This has been exten- 
sively done in Ill-nitrides both in epilayers and quantum 
confined heterostructuresPSl l n this case the phonon- 
assisted transitions "replicate" the ZPL shape and con- 
tain the same information on the exciton localization dy- 
namics. Depending on the localization degree a thermal 
regime can eventually be reached, which requires tens or 
hundreds of picoseconds. 

On the other hand, in samples where exciton localiza- 
tion is weak, the free exciton (FE) thermalization issue 
is relevant and from the study of the phonon- replica emis- 
sion we can expect extracting useful information on the 
exciton distribution over an extended k-range. More- 
over, even if the ZPL shows bound exciton emission (in 
non-intentionally doped (nid) GaN epilayers excitons are 
mainly bound to donor states), the PR bands mostly 
originate from intrinsic excitonic states, making negligi- 
ble the contribution of the extrinsic band. The PR bands 
contain information on the whole (k = and k^O) ex- 
citon population, since the momentum conservation rule 
is assured by the presence of one or more phonons. In 
particular longitudinal optical (LO) phonons, of energy 
ElOi are very effective in determining intense PR-PL 
bands, which are denominated nLO band depending on 
the number n of LO phonons involved in the emission 
process. In fact it has been shown that the shape of the 
2LO PR of the FE band is proportional to the FE popu- 
lation, while the 1LO PR of the FE band reflects the FE 
population times the exciton kinetic energy!^ Therefore 
phonon replicas have been investigated in different semi- 
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conductors to assess both the exciton-phonon interaction 
and the FE dynamics at large k, which is of interest not 
only for fundamental physics but also for applications P 
In particular, the study of the phonon-assisted radiative 
recombinations provides significant information on the 
exciton dynamics in II- VI and III-V nitride semiconduc- 
tors, which show pronounced phonon assisted emission P 
In III-V nitride structures, m ost of the reports on PR 
refer to the steady-state regimeP^We recently showed^ 
that a close comparison with existing models^ for the 
PR emission in GaN epilayers is relevant provided a free 
exciton thermalized distribution is achieved and taking 
into account the complex nature of the exciton band in 
GaN. In this way a full reconstruction of the exciton 
density of states is made possible and the experimen- 
tal PR energy shift, lineshape and intensity are nicely 
reproduced by this model. The k = exciton dynamics 
has been extensively investigated both at the fs and ps 
time scale. Typically a 150 ps radiative recombination 
time and fast dephasing times (equal to 1-2 ps in high 
quality samples at low excitation density) are found for 
the A and B excitonsP^K^ Nevertheless, despite the fact 
that in the 1-100 ps range the exciton dynamics is domi- 
nated by large-k excitons, the PR recombination kinetics 
has been scarcely investigated mainly focussing on the 
quasi-thermal regime.^ But even if the exciton-acoustic 
phonon interaction is quite strong, the thermal regime 
is reached after a finite time. This implies the presence 
of a non-thermal regime, which has not been experimen- 
tally investigated yet. Only theoretical results concerning 
the exciton formation and the approach to thermal equi- 
librium are reported in literature.^^ In this study we 
aim at filling this gap and therefore we explore both the 
non-thermal and the thermal regime. From picosecond 
time-resolved experiments performed on nid GaN epilay- 
ers grown on c-plane sapphire substrates, we show that 
the achievement of a hot thermal distribution of excitons 
requires a few tens of ps, depending on the excitation 
energy, the lattice temperature and the excitation power 
density. The heating process is also studied by resonantly 
exciting the k = FE (Xa exciton). As a consequence 
of the different wavevector intervals explored by the PR 
and the ZPL emission, we provide insights on the exciton 
dynamics in an extended region near the T-point of the 
first Brillouin zone. 



II. EXPERIMENTAL 

The investigated samples are a 11 (im thick nid 
wurtzite GaN epilayer grown by hybrid vapor phase epi- 
taxy (HVPE) on a c-plane sapphire substrate^ (sample 
A) and a 3 /im thick nid wurtzite GaN epilayer grown 
by metal organic vapor phase epitaxy (MOCVD) on a c- 
plane sapphire substrate (sample B). The threading dis- 
location density is as low as 2 x 10 8 cm~ 2 in the HVPE 
sample and lower than 1 x 10 9 cm -2 in the MOVPE epi- 
layer with a not intentional n-type doping w 10 17 /cm~ 3 . 



\ V 




3.2 3.3 3.4 3.5 

E(eV) 

FIG. 1. PL and reflectivity spectrum measured at 15K show- 
ing the details of the ZPL and phonon replica (PRE) bands 



Measurements were performed with the sample placed in 
a closed-cycle cryostat varying the temperature from 10 
to 50 K for time-resolved experiments. Reflectivity mea- 
surements at normal incidence were carried out using a 
CW Xe lamp. Time-integrated photoluminescence (TI- 
PL) measurements were performed under non-resonant 
excitation by a frequency-doubled ps dye laser at 300 nm, 
and the collected light was detected by a cooled silicon 
charge coupled device camera after dispersion through 
a 50 cm flat-field spectrometer with a spectral resolu- 
tion of 1 mcV. Time-resolved photoluminescence (TR- 
PL) measurements were performed using a frequency- 
doubled mode-locked Tksapphire laser pumped by a con- 
tinuous wave Ar + laser providing 1.2 ps pulses at a rep- 
etition rate of 81 MHz. The PL was dispersed through a 
30 cm flat-field monochromator and detected by a streak 
camera apparatus with a 3 ps time resolution. In Fig. [I] 
typical TI-PL spectra taken at T = 15 K for the ZPL and 
the PR region are shown for sample A. Similar spectra 
for sample B are shown in Ref. 'LUJ The main spectral 
features observed in the ZPL region at 10 K are related to 
the neutral donor-bound exciton recombination (D°Xa) 
and to the Xa and X^ exciton radiative recombinations. 
The emission from the 2S state of Xa excitons is also de- 
tected at higher energy. The reflectivity spectrum (blue 
line in Fig. [I]) shows clear excitonic resonances confirm- 
ing the labeling of the PL structures. In this case the 
Xc excitonic resonance is detected at ~ 3.52 eV. In the 
PR region the reflectivity spectrum shows only interfer- 
ence fringes related to the epilayer thickness. From the 
comparison of Fig. [T] and the spectra reported in Ref. 
rnjlit turns out that samples A and B show very similar 
features and pronounced excitonic resonances. The mea- 
sured Stokes shift for the Xa exciton, detected at 15 K, is 
also very similar and is lower than 1.5 meV in both sam- 
ples indicating a weak exciton localization at low T . The 
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difference observed in the exciton energies, as extracted 
from the reflectivity spectra, comes from the different 
residual biaxial strains of the two epilayers and agrees 
with data and calculation already published! 19 * 20 * For a 
detailed discussion on the TI-PL spectra we refer the 
readers to Ref. 10] We only mention that the TI-PL line- 
shape almost coincides with the long delay TR-PL line- 
shape and that it well agrees with the thermal model of 
Segall-MaharP when the discrepancy commonly reported 
for the temperature behavior of the phonon-replicas is 
described by accounting for the complexity of the GaN 
excitonic bandP^ Here we will focus on the time-resolved 
spectra to investigate the typical time scale for reaching 
a quasi-thermal distribution for which the model of Ref. 
[5] can be applied. Since the experimental results turn 
out to be very similar for samples A and B , most of 
the data we will show will refer to sample A. In Fig. [2] 
typical TR-spectra are shown for (a) the ZPL, (b) the 
first LO phonon replica (1LO) and (c) the second LO 
phonon replica (2LO) acquired when exciting the sample 
with 3.55 eV photons (i.e., photons with an excess en- 
ergy of nearly 60 and 35 meV with respect to the bottom 
of the exciton and free carrier band, respectively, which 
is however still lower than the energy £xo = 92 meV 
of the dominant LO phonon). Here, and in the follow- 
ing, the zero time corresponds to the maximum of the 
excitation pulse as measured from the Rayleigh or Ra- 
man scattered light. It appears that for the ZPL most of 
the dynamics involving free exciton states, both A and 
B ones, occurs during the first 40 ps following the exci- 
tation. The modifications seen at later times in the ZPL 
lineshape come from the increasing contribution of the 
donor- and acceptor- bound exciton states. Concerning 
the time evolution of the replica bands we will first focus 
on the 2LO-phonon band whose shape and intensity is 
directly proportional to the exciton population along the 
k-statesPAt early times, the emission is blue-shifted and 
broadened with respect to that measured at delays longer 
than 100 ps. Besides, for delays longer than 40 ps, the 
spectrum gets narrower and shifts to the energy position 
expected for a thermal exciton population. The complex- 
ity of the exciton dynamics, which involves a non-thermal 
and a quasi-thermal regime at the picosecond time scale, 
is suggested by the different time evolutions of the spec- 
tral components of the 2LO-phonon band reported in Fig. 
[3]ja) and Fig. [3](b) , respectively. As a general trend, the 
rise and fall time get shorter with increasing energy of the 
recombining state. The risetime of the high energy states 
(where only FE are present), at the limit of our exper- 
imental time resolution, indicates an exciton-formation 
time «10 ps followed by a PL decay occurring on a 100 
ps time scale, which originates from both the population 
thermalization and decay. At early times the high energy 
side of the phonon replica is faster than the Xa emission 
and we recover a common exponential decay at longer de- 
lays, indicating the establishment of the thermal regime. 
The lower energy states exhibit a slower rise and decay 
as a consequence of the population dynamics and/or the 



exciton cooling. 

In Fig. [3^c) , we show the time evolution of the spectrally 
integrated free exciton emission of the 2LO-phonon band, 
excited with an excess energy of «60 meV for two dif- 
ferent power densities (10 Wcm~ 2 and 100 Wcm -2 ) nor- 
malized to the excitation power density. The fit with an 
exponential rise and decay is reported on the experimen- 
tal data and, for clarity, the experimental time response is 
also shown. In this case the rise and decay times reflect 
the time evolution of the total exciton population free 
from thermalization effects. In the following we will refer 
to the rise-time of the spectrally integrated PL as to the 
" exciton-formation time" . Figure [3] (c) clearly shows a 
decrease in the exciton formation time (and a 20% super- 
linearity) when increasing the excitation power density 
from 10 to 100 Wcm~ 2 . In addition the exciton forma- 
tion time only sligthly depends on the excitation energy 
as supported by the decrease from 13 to 9 ps when re- 
ducing the excitation excess energy below the free carrier 
absorption edge (from 60 to 20 meV with respect to the 
bottom of the exciton band). This suggests, as we will 
discuss later, that we mostly create electron-hole (e-h) 
correlated pairs rather than free carriers, even when we 
excite the samples a few tens of meV above the free car- 
rier continuum. 

Let us now move to the study of the exciton relaxation 
process. In Fig. [4] we show the comparison between 
the 2LO phonon replica time-resolved spectra taken at 
t = 114 ps and t — 10 ps delay time after the excitation 
pulse. It is worth pointing out that the overall exciton 
dynamics is almost independent from the excess energy 
of the exciting photons (not shown here); therefore in 
the following the experimental data will refer to the ex- 
citation energy of 3.55 eV. The exciton relaxation dy- 
namics has been commonly interpreted as a cooling pro- 
cess assuming that an exciton ther mal d istribution can 
be reached within a few picosecond d 14 * 21 * and then that 
we can reproduce the PL lineshape with a Boltzmann 
distribution with a temperature dependent on the delay 
time. In fact, in agreement with previous reports,^ we 
find that, within the experimental sensitivity, the high 
energy tail of the emission is mostly exponential even at 
very short delay times. However, we can perform a more 
accurate analysis by using the fact that the lineshape of 
the 2LO phonon band is directly related to the popula- 
tion along the exciton band by 

I 2 Lo(E,t)K P (e)f(e,t)Q(E-E ) (1) 

where e = (E — E ) is the exciton kinetic energy, E = 
Ex A — 2Elo, p(e) represents the exciton density of states 
and f(e,t) is the exciton distribution function at time t. 
Generally speaking, /(e, t) can even be a non-thermal dis- 
tribution and, in fact, reasonable fits to the time-resolved 
2LO-phonon emission spectra with Boltzmann functions 
for /(e, t) are obtained only at long delays. In Fig. |4j we 
compare the time-resolved 2LO-phonon spectra with the 
corresponding fits using a Boltzmann-like distribution. 
The spectrum at a delay time of 114 ps is nicely repro- 
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duced by p(e) given in Ref. HU]and k B T L = 1.97 meV, T L 
being the temperature of the lattice. The best fit to the 
data at a delay time of 10 ps requires a hotter thermal dis- 
tribution (fcfiT = 3.77 meV). However the quality of the 
fit is quite poor. In particular the states near the peak en- 
ergy reveal a lack of population with respect to the ther- 
mal condition, suggesting a non-thermal exciton distribu- 
tion. In order to quantitatively discuss the exciton ther- 
malization we, first of all, normalize each time-resolved 
spectrum of the 2LO-phonon PL band with its spectrally 
integrated intensity. This operation eliminates the ex- 
citon formation and recombination from the lineshape 
analysis, which, to a first approximation, only affects the 
total number of excitons. The resulting normalized spec- 
tra l2 LO (E,t) is then proportional to the normalized ex- 
citon population N(E, t). The next step is to bypass the 
problem of modelling p(e). To this aim, considering the 
ratio n(E,t) = I% h0 (E,t)/I% h0 (E,t = 200ps) between 
the normalized spectra, we have experimental data that 
do not depend on the actual p(e). In fact, according to 
the experimental data, we can assume a thermal distribu- 
tion at the lattice temperature for t > 100 ps ( hereafter 
t = oo ) . Then from Eq[l] the 2LO-phonon normalized 
emission band at t = oo is given by: 



'2LO 



(E, oo) oc 



p(e)e k B T L 

(fc B T L )(3/2 



e(E-E ) 



(2) 



and therefore from n(E, t) we can extract the exciton dis- 
tribution /(e, t) even if the density of states is unknown. 
In Figs. [5ja) and [5|h) we report the time evolution of 
the normalized spectra l2 LO (E,t) oc N(E,t) at fixed en- 
ergies for two different excitation power densities (10, 
100 Wcm -2 ). The quantity N(E,t) decreases with time 
for large E values, while it increases for small E values. 
These trends are obviously associated with the transfer 
of population towards the bottom of the excitonic band. 
We also observe that the state of energy around 3.310 
eV, nearly corresponding to an excess of energy SksT^/2 
above the bottom of the exciton band, rapidly reaches 



the equilibrium value at temperature T^, i.e., a constant 
value of the population. A faster thermalization is ob- 
served at higher excitation power density, as expected if 
we take into account exciton-exciton interaction. In order 
to get more direct information on the exciton distribution 
f(e,t) we define the quantity j3*{E,t): 



P*(E,t) = ln- 



'2LO 



(E,t) 



z 2LO 



(E,oo) 



In 



(k B T L )Wf(e,t) 



k B T L 



(3) 



j3* (E, t) turns out to be an indicator of the degree of ther- 
malization of the excitons and it is independent on the 
actual density of states. Whenever the hot exciton pic- 
ture applies, i.e., when the exciton distribution f(e,t) is 
a quasi-thermal distribution with a time-dependent tem- 
perature T(t), we get 



P*{E,t) th 



1 



1 



k B T L k B T(t) 



3 / k B T L 
2 \k B T(t) 



(4) 



In the hot exciton limit we therefore expect to get a 
straight line when plotting j3* [E, t) as a function of en- 
ergy with a slope given by \/k B T(t). Any deviation from 
a linear dependence vs. E is therefore a signature for a 
non-thermal regime. In Fig. [6] (a), j3*(E,t) is reported 
for different delay times, in the 2LO-phonon spectral re- 
gion, as extracted from TR spectra. The nominal lattice 
temperature is equal to 10 K and the sample is excited 
with an excess energy of ~ 60 meV with respect to E a 
with a power density of 10 W cm -2 . For clarity we re- 
port the value of j3* (E, t) only in the spectral region of the 
free exciton emission with the low energy side cutoff cor- 
responding to the onset of the 2LO-phonon D°Xa emis- 
sion. In this figure, the energy E is identified with the 
photon emission energy and not only the kinetic energy 
term. We also show the time-integrated PL spectrum of 
the 2LO-phonon replica corresponding to the full thermal 
regime. Strong deviations from a linear trend of /3*(E, t) 
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FIG. 3. PL time evolution at different energies around the 
2LO-phonon Xa exciton (a) in a restricted and (b) over an 
extended time scale, (c) Time evolution of the spectrally in- 
tegrated emission (only the free exciton contribution) of the 
2LO-phonon band compared with the experimental time re- 
sponse (black line) for two different excitation densities. Or- 
ange and red lines are fits as discussed in the text. 



are observed at the early stage of the exciton dynamics 
{t < 40 ps) . In particular a decrease in the exciton popu- 
lation for small kinetic energies with respect to a thermal 
distribution is found. In addition, we also observe a flat- 
tening of /3* {E, t) on the high energy side (corresponding 
to the tail of the PL band) , which could be a further indi- 
cation for a non-thermal distribution. However, this part 
of the spectra is affected by a relevant background whose 
subtraction is delicate making reliable conclusions rather 
hazardous. It turns out that the distribution function 
behaves like a Boltzmann function only at longer delays 
(t > 50 ps) and f3*(E,t) shows a linear dependence over 
the whole energy range covered by the 2LO-phonon band 
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FIG. 4. TR spectra of the 2LO-phonon band at different 
delays along with a fit as discussed in the text. 

with a decreasing slope in agreement with a cooling of the 
exciton population. Finally j3* (E, t) becomes constant 
and vanishes, as expected, when the distribution reaches 
the lattice temperature. Referring to Fig. RsFa), the de- 
viation of /3* (E, t) from the linear behavior observed for 
t < 40 ps indicates the lack of excitonic population not 
only with respect to the thermal regime at T = Tl but 
also with respect to a hot exciton distribution. There- 
fore our findings clearly evidence the presence of a non- 
thermal distribution of the exciton population at short 
delay times. It comes out that the cooling picture of a 
hot exciton distribution does not apply if we refer to the 
whole excitonic band. Moreover, the lack of population 
at the Xa energy with respect to the hot thermal distri- 
bution cannot be ascribed to a localization mechanism. 
In fact, the absence of any significant localization of the 
excitons is proven by the comparison of reflectivity and 
PL spectra for the ZPL (Fig. [IJ where it turns out that 
the Stokes shift at 10 K is 1 meV, a value not large 
enough to explain the experimental results displayed in 
Fig. |6ja). In Fig. |6jb), similar spectra are reported for 
an excitation power density of 100 Wcm~ 2 . No major 
change is found except that the quasi-thermal regime is 
reached for slightly shorter times with respect to the low 
excitation condition, as already shown in Fig. [5] Up to 
50 K the non-thermal regime can still be detected during 
the first 30 ps and similar results for f3*(E, t) are observed 
when we change the photon excitation energy for an ex- 
cess energy lying in the 6 to 60 meV range. A specular 
behavior is instead observed when we resonantly excite 
the k=0 Xa exciton and we study the heating of the cor- 
responding exciton distribution. In Fig. [7]Ja) TR spectra 
are shown for different delay times after the resonant ex- 
citation at E = 3.303 eV slightly below the Xa exciton. 
Apart from the initial Raman signal detected for the first 
10 ps, that dominates the emission spectra at the early 
stage, the exciton line and also that of the D°Xa emission 
are visible at t > 20 ps. In Fig. 0b), we report j3*(E,t) 
for t > 20 ps when the Raman signal can be neglected 
and we can discuss the evolution of the exciton distribu- 
tion without spurious emission. Due to this cut occurring 
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FIG. 5. Time evolution of the normalized exciton population 
N(E) at different energies (equi-spaced) in the 2LO band for 
f0 Worn" 2 (a) and 100 Wcm -2 (b) excitation power densi- 
ties. The continuous (dashed) lines refer to states lying above 
(below) 3.307 eV. 

at the early stage in the time evolution of the excitons 
it is not possible to have a definitive proof of their non- 
thermal distribution under this excitation condition. Still 
P* (E, t) shows a slope opposite to that reported for non- 
resonant excitation, denoting a cold exciton distribution. 
During the thermalization process we observe an exci- 
ton heating with a population transfer between low and 
high energy states until the equilibrium distribution at 
temperature is reached. Quite remarkably we find 
that the time scale for reaching the equilibrium distri- 
bution is almost the same for the cooling and heating of 
the exciton population. Finally it is worth noting that, 
when comparing the PL spectra shown in Figs. [6] and [7| 
a saturation of the acceptor emission is observed (Fig. 
when increasing the excitation power density while it is 
its absence which is observed when exciting the sample 
resonantly (Fig. [7]). 

III. DISCUSSION 

The analysis of the 2LO PR-PL lineshape, without 
any supplementary modeling of the exciton density of 
states, already gives important information on the ex- 
citonic distribution function. However, to quantitatively 
study the thermalization process, we have to analyze and 
compare the time evolution of the mean energy < e(t) > 
and that of the effective temperature T(t) of the exci- 
ton distribution. < e(t) > can be directly obtained from 
the experimental data by considering the ratio between 
the spectrally integrated lLO-phonon emission and the 
2LO-phonon band p2l while T(t), according to the usual 
approach,^ is given by an exponential fit of the high en- 
ergy tail of the 2LO-phonon PL spectra. Obviously for 
a thermal exciton distribution and a standard excitonic 
density of states we expect to have < e(t) >= 3k B T(t)/2. 
The two data sets, < e(t) > and 3fcsT(t)/2, are reported 
in Figs, and[8^b) for two excitation power densities. 
We observe that in the early stage < e(t) > is larger than 
3kgT(t)/2i thus denoting a non-thermal distribution of 
the excitons. At long delay time < e{t) > and 3/csT(i)/2 
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FIG. 6. Time evolution of the quantity f3*(e,t) for different 
delay times for non-resonant excitation: 10 Wcm -2 (a) and 
100 Wcm" 2 (b). The solid black line in (a) and (b) is the 
time-integrated PL spectrum. 
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FIG. 7. (a) Time resolved spectra taken at different delay 
times for resonant excitation, (b) Time evolution of the quan- 
tity /3*(e, t) for the resonant case. For clarity the 2LO-phonon 
PL spectrum in the fully thermalized regime is also shown 
(black solid line). 



reach the same asymptotic values, as expected for an 
excitonic thermal distribution. In Fig. [8jc) the value of 
< e(t) > is reported for two different lattice temperatures 
T L = 10 and 50 K. At T L = 50 K we see that the exciton 
cooling is significantly faster than that at Tl = 10 K. In 
Figs. [8]ja)-[8]jc) the continuous lines are best fits to the 
experimental data obtained with the expression for the 
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FIG. 8. Time evolution of the exciton average kinetic energy 
< e > and the thermal energy along with the fits as dis- 
cussed in the text for the non-resonant excitation case, (a) 10 
Wcm~ 2 . (b) 100 Wcm -2 .(c) 10 K and 50 K, nominal lattice 
temperature at 10 Wcm -2 . 



time-dependent mean energy < e(t) > (or 3kgT(t)/2) 



< e(t) >=< eoo > 



Ae t / T - 1 



with 



> 



A = V<^> + v< 

V< £ > - V< e oo > 



(5) 



(6) 



where < > and < eo > are the equilibrium and ini- 
tial (i.e. at t = 0) mean energy . The thermalization 
time constant r is determined by the exciton-aco ustic 
phonon interaction due to the deformation potentia l 14 * 22 ! 
and does not depend on the lattice temperature Tl . This 



expression has been very often derived assuming a ther- 
malized distribution, therefore referring to the temper- 
ature. Nevertheless it holds independently from this 
assumption]!" The values determined for the time con- 
stant of < e > and 3ksT(t)/2 are respectively r = 23 ±2 
ps and 32±2 ps at 10 Wcm' 2 and t = 15±1 ps and 24±1 
ps at 100 Wcm" 2 ; for T L = 50 K we get r = 12±1 ps. As 
expected the cooling process becomes faster when rising 
either the excitation power density or the temperature. 
In addition the data fitting points out a strong discrep- 
ancy between < e > and 3kgT(t)/2, strenghtening the 
non-thermal distribution character of the excitonic gas. 
We can now compare the energy relaxation time with 
the values predicted by theories. From r = 23 ps (32 
ps), by means of the expression reported in Ref. 1221 
we get an effective deformation potential (DP) constant 
Ei = 20 ± 1 eV (17± 1 eV) to be compared with the value 
Ei = 13±1 eV reported in Ref. [Hand Ei = 9±2eV usu- 
ally found both theoretically and experimentally (e.g. see 
Ref. [53J). Provided the value we find is nearly twice the 
" standard" one, and that the latter would lead to a ther- 
malization time of about 100 ps, we should consider other 
scattering mechanisms in addition to that of the defor- 
mation potential. A contribution to the exciton thermal- 
ization from exciton-exciton scattering is evidentiated by 
the acceleration of the dynamics when rising the excita- 



tion power density from 10 to 100 Wcm 



the exci- 



ton density from rj 10 cm~ J to 
the role of the inelastic exciton-exciton collision in the 
population re distrib ution in GaN has been reported by 
several autorsP^EH An evaluation of the exciton-exciton 
scattering rate can be obtained by the dephasing time. 
Indeed, according to the results of Pau et al. (see Ref. 
|21)|) this contribution to the lincwidth depends on the ex- 
citon density Nx as Txx = Pxx^bEb^x , where &b is 
the exciton Bohr radius, Eb the exciton binding energy 
and /3xx a dimensionless efficiency factor which for GaN 
is between 5 and 7. With &b — 3 nm and Eb = 24 
meV, as found from an hydrogenoid model, we get a de- 
phasing time T2 = irTxx/h in the range 20 to 40 ps for 
Nx = 10 16 cm -3 . We get instead a value of 2 to 4 ps 
for Nx = 10 17 cm -3 , which is an order of magnitude too 
short to account for the experimental results. From our 
data, it turns out that, in the thermalization process, the 
inelastic exciton-exciton scattering contributes but with 
an efficiency ten times lower than in the dephasing pro- 
cess, where collisions are prevailing without any change 
in the energy of the single exciton. It is also worth point- 
ing out that the exciton-free carrier scattering process 
can in principle also contribute to the exciton relaxation 
but, from time-resolved spectra, the number of free carri- 
ers seems negligible with respect to that of the excitons. 
Therefore other mechanisms have to be considered and 
in particular the piezoelectric interaction (PZ) with the 
acoustic phonons as well as the impurity (defects) inelas- 
tic (D°Xa) scattering. The role of the PZ interaction in 
the exciton thermalization process in GaN epilayers has 
been extensively investigated by Kokolakis et al. (see 
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Ref. [T5))and it has been found nearly equal to that of the 
DP interaction for low residual doping. Due to screen- 
ing effects it becomes negligible for a 10 17 cm -3 impurity 
content. Provided that the PZ mechanism does not give 
a good fit to the experimental data we consider its con- 
tribution to the thermalization as being marginal. When 
considering the scattering from neutral don ors th e im- 
pact rate is roughly given by Tdx = a r v gND I2ZI3S1 where 
Nd is the neutral donor concentration and v g = ^/2e/m 
is exciton group velocity, being e and m the exciton ki- 
netic energy and mass. Then, for example, with e = 2 
meV and Nd = 10 17 cm~ 3 we get a very short collision 
time tdx = (Fdx) -1 ~ 400 fs. Therefore, even if most 
of the collisions are of elastic nature and that they only 
contribute to the exciton dephasing process,^ we can ex- 
pect a relevant contribution to the exciton thermalization 
from the inelastic processes as already found in GaAs 2 ^ 
and InPpS This mechanism is also in agreement with the 
observed behavior when changing the lattice temperature 
Tl provided the scattering rate only depends on the ex- 
citon energy and not on T^. Thus, as observed in Fig. 
|8jc), the initial energy relaxation towards the exciton 
thermalization is independent of Tl. Since the exciton 
initial excess energy is identical for T = 10 and 50 K, 
the acceleration observed in Fig. [8^c) is only becoming 
apparent as a consequence of the higher final energy of 
the excitons at 50 K. This does not hold for the DP (PZ) 
interaction. In fact the fit of the experimental data of 
Fig. [8|c) with Eq. [5] requires a decrease in the time con- 
stant related to the DP (PZ) interaction, when rising the 
temperature, in disa greem ent with its expression, which 
is independent of TzJ 14 ' 22 1 

The previous discussion deals with the time evolution of 
a hot exciton population that has reached a finite and 
constant value and evolves towards the thermal equi- 
librium with the lattice. Let us now briefly comment 
on the early stage of the exciton dynamics when the 
states near the bottom of the excitonic band become 
populated, i.e., the exciton formation stage, which in our 
measurements corresponds to the time interval between 
the arrival of the excitation pulse and the observation 
of the maximum of the spectrally integrated PL of the 
2LO-phonon replica(Fig. [3](c)). The corresponding PL 
rise time turns out to be ~ 10 ps, nearly independent 
on the energy of the excitation photon, which is only 
slightly shortened when rising the excitation power den- 
sity. When exciting the sample « 60 meV above the 
state, i.e 35 meV above the free carrier continuum, we 
create both free carriers and correlated e-h pairs. Free 
carriers relax in energy and form excitons via a bimolcc- 
ular (BM) mechanism with a coefficient B = 1.2 x 10~ 18 
ps -1 cm 3 PSl Correlated e-h pairs reach the bottom of 
the exciton band losing their excess energy by means 
of inelastic interactions with acoustic phonons, impuri- 
ties and defects. At the excitation power densities used 
in our experiment the bimolecular generation only con- 
tributes at the highest power density (100 Wcm -2 ), i.e. 



10 17 carrier/cm 3 ), providing an exciton formation time of 
12 ps, not far away from the observed rise time of the 
2LO-phonon replica. In addition the slight supcrlinearity 
shown in Fig. |3jc) well agrees with the BM mechanism. 
In our experiment at 10 Wcm" 2 the exciton formation 
time is only sligthly longer and is still in the range of 
10 ps. Then it turns out to be shorter than the time 
needed to reach the bottom of the exciton band when an 
energy loss of 30-40 meV by emitting acoustic phonons is 
considered. Moreover the BM generation would provide 
a risetime of « 100 ps.^ Therefore, also for the exciton 
formation, inelastic collisions with neutral donors appear 
as the most likely process for a fast exciton formation, or 
more precisely, for a fast relaxation of correlated e-h pairs 
towards the bottom of the exciton band. One inelastic 
collision, giving rise to the donor ionization, is sufficient 
to produce a 20-30 meV loss in the exciton energy, thus 
providing a fast transfer from the high to the low energy 
exciton states. The fast rise in the ZPL signal, equal to 
that reported for the 2LO-phonon PL in Fig. [2] at low ex- 
citation power density is mainly due to exciton-impurity 
scattering, which avoids the bottleneck resulting from the 
reduced efficiency for the exciton-acoustic phonon scat- 
tering. 

IV. CONCLUSIONS 

From a careful analysis of low temperature time- 
resolved spectra of the phonon replica emission in GaN 
epilayers we have evidenced the transition from the non- 
thermal to the thermal regime of the exciton population. 
Our experimental results provide clear signatures that 
the processes leading to the exciton formation and ther- 
malization are relatively slow when the exciton density 
is lower than 10 17 cm" 3 . The establishment of a thermal 
Maxwcll-Boltzman distribution for the exciton popula- 
tion is observed after sa 50 ps, mainly as a consequence 
of inelastic collisions with impurities rather than with 
acoustic phonons. From the comparison with models and 
predictions for the deformation potential, piezoelectric 
interaction and impurity scattering we show that neu- 
tral donor scattering can account for the exciton forma- 
tion and relaxation. Only thanks to the residual doping, 
the bottleneck regime for the radiative excitons at k = 
is avoided. Nevertheless, for an exciton density > 10 17 
cm~ 3 , the exciton-exciton collisions become relevant in 
the exciton formation and thermalization, contributing 
to the cooling via inelastic processes. 
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